Duality and integer quantum Hall effect in isotropic 3D crystals 
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We show here a series of energy gaps as in Hofstadter's butterfly, which have been shown to exist 
by Koshino et al [Phys. Rev. Lett. 86, 1062 (2001)] for anisotropic three-dimensional (3D) periodic 
systems in magnetic fields B, also arise in the isotropic case unless B points in high-symmetry 
directions. Accompanying integer quantum Hall conductivities (o- xy ,a yz ,a zx ) can, surprisingly, 
take values oc (1, 0, 0), (0, 1, 0), (0, 0, 1) even for a fixed direction of B unlike in the anisotropic case. 
We can intuitively explain the high-magnetic field spectra and the 3D QHE in terms of quantum 
mechanical hopping by introducing a "duality", which connects the 3D system in a strong B with 
another problem in a weak magnetic field (oc 1/B). 
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I. INTRODUCTION 

Appearance of energy gaps in three-dimensional(3D) 
systems is a hallmark of quantum mechanics in crystals. 
There, a periodicity, be it atomic or density-wave for- 
mation from some mechanism, gives rise to a Bragg's 
reflection accompanied by energy gaps as Bloch's theo- 
rem dictates. Apart from this, there is few occurrences 
of energy gaps in 3D (unless of course non-perturbative, 
many-body states such as the BCS state are involved). 

In two-dimensional (2D) systems, by contrast, we do 
have a remarkable realization of gaps when a magnetic 
field is applied. For a free 2D electron gas the spec- 
trum even coalesces into a series of sharp Landau levels. 
When the 2D system has some periodicity (thought of 
as arising from a periodic potential or tight-binding ar- 
ray of atoms), application of a magnetic field gives rise 
to a fractal energy gaps, which is called "Hofstadter's 
butterfly" 1 . One way to realize why the butterfly ap- 
pears is that the two quantizations, one being Bloch's 
band structure and the other the Landau quantization, 
interfere with each other. The interference should occur 
when the band gap, ft 2 /(ma 2 ) (to: electron mass, a: lat- 
tice constant), and the cyclotron energy, hu) c , are similar, 
which indeed amounts to the condition for the appear- 
ance of the butterfly, <f> ~ </>o , where <f> is the magnetic flux 
penetrating a unit cell and 4>o = h/e the flux quantum. 

Once we have an energy gap in a 2D system it is known 
that the integer quantum Hall effect (IQHE) should oc- 
cur. Thouless et aP in fact have derived a general for- 
mula for the quantized Hall conductance in 2D periodic 
systems without any assumption on the periodicity, and 
shown that the Hall conductance can be written in terms 
of some topological invariants. 

A natural question is whether we can extend such ar- 
guments to 3D systems. It has been shown that, if 
an energy gap exists in 3D, it should accompany an 
IQHE where each component of the Hall tensor er^ is 
quantized 3-6 . A big question, then, is whether and how 
gaps can appear in 3D. We have previously shown that 
the butterfly-like spectrum indeed appears in anisotropic 
3D (quasi- ID) lattice and have derived the quantized 
Hall tensor for each gap. 12 For the anisotropic case we 



can intuitively understand these in terms of a mapping 
between the 3D case and a 2D case. For the isotropic 
case, however, it seems difficult to have the energy gaps 
since the magnetic subbands are likely to be overlapped 
by the large dispersion along the magnetic field. Sev- 
eral authors 7 ' 8 have computed the energy spectra of the 
isotropic 3D tight-binding system under magnetic fields 
pointing to high-symmetry crystallographic directions, 
and Kunszt 8 concluded that there is in general no en- 
ergy gaps although the possibility of gaps cannot be ex- 
cluded. So the existence of gaps and IQHE remains an 
open question. 

In this paper, we show that we do have energy gaps 
in isotropic 3D lattices rather universally in that all we 
have to have is the sufficiently large magnetic field B 
pointing to general directions (i.e., off the high-symmetry 
axes). We first show this by calculating energy spec- 
tra for various directions of B. We have also calculated 
the Hall tensors, o~ xy , a yz , a zx , when the Fermi energy is 
in each energy gap. They turn out to have values, e.g. 
(1, 0, 0), (0, 1, 0), (0, 0, 1) in an appropriate unit, even for 
a fixed direction of B, which is surprising and unlike in 
the anisotropic case. 

We then intuitively explain such structures in the spec- 
tra by introducing a "duality" , which connects the 3D 
system in a strong B with another problem in a weak 
magnetic field (~ 1/B). This enable us to regard the 
band structure and the quantization of the 3D Hall con- 
ductivity as graphically arising from quantum mechani- 
cal hopping between "semiclassical" orbits in the weak-f? 
case. 



II. FORMULATION OF 3D BLOCH ELECTRONS 
IN MAGNETIC FIELDS 



We take a non-interacting tight-binding model in 
a uniform magnetic field B — (B X ,B V ,B Z ) pointing 
to an arbitrary direction. Schrodinger's equation is 
~^2ijtij e i i } j = where ipi is the wave function 

at site i, the summation is over nearest-neighbor sites, 
and 9ij = j- J* A ■ dl is the Peierls phase with A the vec- 
tor potential, V x A = B. We consider a simple-cubic 
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lattice with a lattice constant a. Following Kunszt and 
Zee, 8 we take a gauge A = (0,B z x — (B Z B X /B y )(y — 
a/2),B x y — B y x), where z is cyclic. We can then write 
ipimn — e lX ' n F^ m , where l,m,n are site indices along 
x,y,z, respectively. Schrodinger's equation is then re- 
duced to a two-dimensional tight-binding model, 

-t Fi 1 — t -t e 2,Ki <t>'\!—(<t>x/<t>y){m--lj\ p 

t'X- 1 t — l,m v x -L £+l,m h y^ ± L : m — 1 

_. -2ni(t> z [l-(4> x /<t> y )m]p 

- 2t z cos[27r(0 x i - 4> y m) + X z }F Wm = EF lim , (1) 

where t; L {i — x,y,z) is the transfer energy and = 
Bi<x /4>q with (f)Q = h/e is respective number of flux 
quanta penetrating the facet of the unit cell normal to 
ej. Here we consider rational fluxes, 

{<j>x,4> y ,<Pz) = $ x (n x ,n y ,n z ), $ = P/Q (2) 

with P, Q: integers and n,'s are mutually prime integers. 
By introducing j = n y l — n x m, we can eliminate m from 
the phase factors (exponentials and the argument of co- 
sine in Eq. (1)), so we can write Fj, m = e tXyTn Gj, where 
G is determined from a one-dimensional equation, 

txGj+ ny txGj— Uy 

+ „ i [2'r*-^-0'+nx)-A H ] r , . -i(27T<S> ^- j - \ y ) „ 

bye y ^j+n^ vy" ^j — nx 

- 2t z cos(-27r$j - \ z )Gj = EGj. (3) 

Since the phase factors have a common periodicity 
n y Q, we can apply the Bloch-Floquet theorem to have 
Gj+nyQ = e lXxrlyQ Gj, and the Hamiltonian is reduced to 
a n y Q x n y Q matrix. 



III. ENERGY SPECTRA AND HALL 
CONDUCTIVITIES 

Our aim is to search systematically for gapful spectra 
for the isotropic case with t x = t y = t z = t and to de- 
termine the Hall conductivities for each energy gap. We 
have numerically solved Eq.(3) and obtained the energy 
spectra versus Fig. 1 shows the results for typical field 
directions: (n x ,n y ,n z ) — (1, 2, 3), (1, 1, 2), (1, 1, 1), and 
(0, 1, 2), which differ in symmetry. In each case we took 
16 points for each of the Bloch wavenumbers X x ,X y ,X z . 

We can immediately see that a series of gaps appear for 
(n x ,n y ,n z ) = (1,2,3), while otherwise we have at most 
solitary gaps. From such results we have found that the 
spectrum has a series of gaps when no two rii 's coincide 
and Hi ^ 0. So we can say that the gapful spectrum can 
be expected when B points to general directions (i.e., off 
high-symmetry crystallographic axes). The gaps shrinks 
when B is too small (Ba 2 /<po -C 1) obviously. 

Now we proceed to the quantum Hall effect. As men- 
tioned, Ep in an energy gap dictates 3-5 quantized Hall 
conductivities, 



{o- yz ,<j zx ,<j xy ) = - — (m x ,m y ,m z ), (4) 

where raj's are integers satisfying a Diophantine's equa- 
tion, 

v B = s + m x 4> x +m y (j)y+m z (j) z . (5) 

Here s is an integer, vb is the filling of the tight-binding 
band, and raj's assigned to each gap are topological in- 
variants, i.e., they never change when external parame- 
ters (magnetic field, transfer energies, etc) are changed 
as long as the gap remains. So Eq.(5) tells us that, if 
we consider vb as a function of (4> x , 4> y , <f> z ) in a 4D pa- 
rameter space (i>b , 4> x ,4> y ,(j) z ), the v& plotted along a gap 
assigned with (m x ,m y , m z ) will form a '3D plane' having 
a gradient (m Xl m yi m z ) in the 4D space. 

In fact we can readily translate the E-$> diagram in 
Fig. 1 to the vb-Q diagram by counting the number 
of subbands below each value of E. There, the gradi- 
ent dvB/d® for each gap gives m x n x + m y n y + m z n z . 
We can then scan the direction of B to determine all 
of m x , m y , m z by keeping track of each gap in ques- 
tion, which is exactly how we have obtained the Hall 
integers here. The result is shown as triple integers in 
Fig. 1. We can see that the largest series of gaps have 
(m x ,m y ,m z ) = (±1, 0, 0), (0, ±1, 0) or (0,0, ±1), where 
only one of the Hall components is nonzero. This is strik- 
ingly new in that one and the same butterfly contains all 
the cases with m x ^ 0, m y ^ 0, m z ^ 0, respectively, for 
a fixed direction of B. This is in sharp contrast with the 
butterfly in the anisotropic case 12 where one component 
a yz is identically zero with x being the most conductive 
direction. Hence the versatile behavior of (m x ,m y ,m z ) 
may be regarded as a hallmark of the isotropic 3D QHE. 

When the above condition is not satisfied (i.e., when 
two nj's coincide as in (1,1,2) and (1,1,1) or at least 
one m = as in (0, 1, 2)), only solitary gaps or zero gaps 
(marked with dashed lines) appear. Hasegawa has shown 
for the case of (1, 1, 1) that two bands touch where the 
energy dispersion is similar to those of a Dirac particle. 
We have similar band touching in Fig.l(b)(l, 1, 2) as well, 
for which a detailed calculation shows that the dispersion 
around the band touching is indeed linear in X x , X y , X z . 

The reason for the band touching can be accounted 
for as follows. For (n x ,n y ,n z ) = (1,1,2) (b), the gap 
for (m X7 m y ,m z ) — (m, 0,0) and that for (0, ra, 0) ap- 
pear along the same line in the vb-$ diagram, since the 
problem is symmetric against x <-> y. If there were 
no band touching, the Hall integers for (ra, 0, 0) and 
(0, m, 0) would become indefinite, which contradicts with 
the topological argument that we can always determine 
raj's uniquely when the gap is nonzero. So we can con- 
clude that the band touching always occurs when two 
rij's coincide. This explains why the series of gaps seen 
in case (a) disappear as we make n x — n y in (b) and 
then n x = n y = n z in (c). We still have to explain why 
the butterfly does not appear in case (d) with one of n^s 
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being zero. We address this question from a different 
viewpoint in the next section. 



3 



(a) (n x ,n y , n z ) = (1,2,3) 
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(b) (n Xl n y , n z ) = (1.1,2) 
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FIG. 1. Energy spectra for the isotropic simple-cubic lat- 
tice in the magnetic flux ((f) x , <j> y , 4> z ) — $ x (n x ,n y ,n z ) (as 
depicted in the inset). Each spectrum is plotted against $ 
with the direction of B fixed at (a,)(n x ,n y ,n z ) = (1,2,3), 
(b)(l,l,2), (c)(l,l,l), (d)(0,l,2). Triple integers represent the 
Hall integers, (cr xy , <r yz , u zx ) in units of (—e 2 /ha). Dashed 
lines in (b),(c) delineate the zero-gap. 
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IV. DUALITY FOR B<->l/B 

Electrons in magnetic fields arc usually described in 
fc-space. A textbook example is the semiclassical pic- 
ture, where an electron is treated as a classical parti- 
cle driven in fc-space by a weak magnetic field. As B 
becomes stronger, the semiclassical orbits begin to be 
mixed and the picture breaks down. Actually this quan- 
tum mixing plays a critical role in opening the gaps and 
in quantization of the Hall conductivity as shown in the 
following. 

We propose here that there exists a duality that relates 
a semiclassical picture in fc-space (weak-B case) to an- 
other semiclassical picture in the opposite limit of strong 
magnetic field, in which the guiding center of the cy- 
clotron motion drifts along equipotential contours in the 
real space. Using a quantum mechanical mapping from 
the weak-i? case to the strong- £? case (where B translates 
into 1/B), we are able to discuss the mixing of semiclas- 
sical orbits in our problem (3D lattice in B) in a clear 
cut. 

Let us derive the duality in 2D and then extend it to 
3D. We consider a 2D Bloch electron with the dispersion 
£{k x , ky) in a magnetic field B perpendicular to the plane. 
If we define the dynamical wavevector K = k + (e/h)A, 
the Hamiltonian in /c-space is written as 

H = e(K x ,K y ), 

[K x ,K v ] = -ieB/n. (6) 

This is an exact quantum equation except that we have 
neglected the inter-band mixing due to the magnetic 
field, which is allowed when the periodic potential is 
strong enough (or the magnetic field is weak enough). 
The equation reduces to the semiclassical one in the limit 
eB/h -> 0. 

If we consider another 2D system in a magnetic field 
B with a periodic potential V(r) with r = (x,y), the 
Hamiltonian in real space is written as 

H = ^-(p + eAf + V(r) 

4 n2 + y (^^ Xl1 )' (7) 

where II = p + eA and R = r — £ is the cyclotron- 
motion guiding center with the relative coordinate given 
as £ = -(l/eB)e z x II. 

When the magnetic field is strong enough (or the po- 
tential V is weak enough) so that (i) the magnetic length 
I = y/K/{eB) is much shorter than the length scale over 
which V varies, and that (ii) the different Landau levels 
are not mixed, we can consider II to be frozen. More 
precisely, the kinetic energy (the first term in Eq.(7)) re- 
duces to a constant from (ii), while the second term can 
be approximated as V(R) from (i). The Hamiltonian 
then becomes, up to a constant, 



HnV(X,Y), 

[X,Y]=ih/(eB). (8) 

We can see that Eqs. (6) and (8) are identical if wc 
translate as 

K x *» v > T~7E- (9) 

So we can see that there exists a duality between the 
2D system with a strong periodic potential in a weak B 
(one band tight-binding model) and the 2D system with 
a weak periodic potential in a strong field (one Landau 
level approximation). For the Landau level in the latter 
case we expect that the duality holds best in the lowest 
Landau level, since the wavefunction is most localized in 
that level to make Eq. (8) valid. 

The duality introduced here provides a physical basis 
for Hofstadter's observation 1 on the mapping between 
the tight-binding and weak-potential systems 9 for the 
square lattice, and also generally applicable to any 2D 
lattices, or, remarkably, 3D cases as well as we shall 
see below. Another duality is found by Ishikawa et a/ 13 , 
which is distinct from ours in that one Landau level is 
mapped to one Landau level in the latter while one Lan- 
dau level is mapped to one Bloch band in the present 
case. 

Now we apply the duality to 3D systems, which en- 
ables us to have a fresh look at the 3D QHE. We take 
the Bloch electron with the dispersion e(k x , k y , k z ) in a 
magnetic field B = (B x ,B y ,B z ). Let us again denote 
the direction of B as n = (n x ,n y , n z ) whose components 
will be assigned a set of integers. We set a new frame 
(x', y', z') in which the magnetic field point to (0, 0, 1). If 
we take the vector potential A with A z > = 0, we have 

H = e(K x r,K y >, k z >), (10) 

where Ki = fcj + (e/K)Ai (i = x',y') with [K x i,K y i] = 
—ieB/h while k z > remains a quantum number. We 
can thus apply the above duality to this system, which 
implies that the 3D system can be mapped to a sys- 
tem in a fictitious 2D space (K x i , K y i ) with a poten- 
tial e(K x i, K y i, k z i) in a fictitious magnetic field B* with 
U/(eB*) = eB/h. The mapping is valid for a large B* , 
i.e., a small B. The fictitious potential is a cross section 
of the 3D dispersion incised at k z > by a plane perpen- 
dicular to B. As we shift k z i, the 2D potential changes 
accordingly, and the spectrum is determined by taking 
all the possible values of k z > . 
We take a cubic dispersion, 

e(k x ,k y , k z ) = —t x cos k x a — t y cos k y a — t z cos k z a, (11) 

which is written in the new frame as, 

e(k x ',ky',k z i) = - t x cos[G x • fc^ + (n x /n)k z >a\ 

— t y cos[Gj, • fc^ + (n y /n)k z ra] 

- t z cos[G z ■ fcj_ + (n z /n)k z 'a], (12) 
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where k± = (k x < , k y '), n = yjn% + n y + n 2 z and 

G x = a(n z n x /n, -n y ), 
G y = a(n y n z /n,n x ), 

G z = d(-(nl + n 2 y )/n,0), (13) 

with a = a/ \Jn 2 + n y . Now we can see that our 3D 

problem is reduced to a fictitious 2D system having three 
periods [while the Hofstadter problem (2D Bloch system 
in B) is reduced to a ID system having two periods]. The 
Hamiltonian for the fictitious 2D system is 

H = ±-(p* + eA*) 2 + e{K x ,,K y ,,k z ,), (14) 

where p* = (d/dK x i, d/dK y >) and A* is a vector poten- 
tial for the fictitious magnetic field (0,0, B*). Since we 
assume a large B*, Schrodinger's equation can be written 
within the basis for the lowest Landau level, which reads, 
after a certain phase transformation, 

-i x * J+nv - t x * 3 -n y -i y e i[2 ^ U+ ^ ) - Xv] iSf j+nx 

- 2t z cos(-2tt<F/ - A z )*j = EVj. (15) 
Here the hopping parameters are 

t i = t i e- G '* im ''\ (16) 

I* = \Jhj (eB*) is the fictitious magnetic length, and 
X y , X z are given by 

X y = -<f> — a H -k z >a 

X z = — k z >a, (17) 

n 

where a labels the wavefunction in the lowest Landau 
level. Now the duality implies that Eq.(15) should reduce 
to Eq.(3) in the strong B* limit (i.e., weak B limit), since 
I* should become smaller than the potential range l/Gi 
so that ti tends to i,. 

V. INTERPRETATION OF THE 3D QHE FROM 
THE DUALITY 

The mapping via the duality enables us to see how 
the total spectrum comes from dispersion relations ver- 
sus k z >, which helps intuitively understanding how the 
energy gaps open and how the Hall conductivities are 
quantized. Since the mapping is limited to the weak-i? 
case, we examine the situation where the gaps begin to 
open in the quantum regime approached from the semi- 
classical one. We show in Fig. 2 the band structure versus 
k z > for the fictitious system (Eq.(15)), along with the fic- 
titious potential e(K x >, K y >, k z > ) for fixed values of k z >. 



Let us first look at how the energy spectra can have gaps 
for general directions of n[oc (1,2,3) in Fig. 2(a)]. e ac- 
commodates closed orbits around its peaks and dips in 
the "semiclassical picture" in the language for the weak- 
B (strong-!?*) case, where the area enclosed by each orbit 
must be a multiple of 2-k/I 2 . So the different wells have 
different sets of discrete levels, where each level moves 
on the energy axis as k z i is changed, since e changes its 
form with k z > ■ When the levels belonging to the wells ad- 
jacent in K x iK y i -space coincide, the states will strongly 
resonate quantum-mechanically and an energy gap will 
arise. In Fig. 2(a), we show the case where the gaps be- 
gin to open with (4> x , 4> y , <fi z ) — (1/9) (1, 2, 3), for which 
the strong mixing of orbits is seen to result in significant 
level repulsions. We attach for comparison the result 
for an almost semiclassical case when the magnetic field, 
(<f) x , 4> y , <j> z ) — (1/45) ( 1 , 2, 3), is so small that the mixing 
is weak and the level repulsion is almost negligible except 
for the middle of the total band. 

In terms of the semi-classical orbits, the change in k z < 
causes the orbits to hop to adjacent positions at every 
resonance as shown in Fig. 2. A virtue of this picture is 
that the distance and the direction of the hopping exactly 
indicates the Hall integers, as shown in the following, If 
we apply an infinitesimal electric field E to the system, 
k z i is dragged adiabatically according to an equation, 
h(dk z i /dt) = —eE z >, where E z > is the component parallel 
to z'. After k z i is changed by Sk z < = 2-K/{an) (the period 
with which the spectrum repeats itself), every state must 
come back to the equivalent position in the reciprocal 
unit cell. Note here that the cell boundaries (represented 
in Fig. 2 as white lines) also move as k z > is changed, since 
the Brillouin zone boundaries are oblique with z' || B. 
The increment in /e-space, SK = (SK X > , 5K y i , Sk z i ) , with 
which the orbit is shifted over the one period satisfies 

SK = =F— dm, 

a 

Sm = Sm x e x + 5m y e y + Sm z e z , (18) 

where — and + correspond to positive and negative E z i , 
respectively, and <5m/s are integers assigned to each sub- 
band. The integers satisfy a relation 8m ■ n = 1, since 
8k z i = 2w/(an). We can immediately translate SK into 
the motion in the real space normal to z' using the re- 
lationship between the relative coordinate £ and the dy- 
namical wavenumber K, 

St = -^e z , x SK. (19) 

Since the system has no dissipation as long as Ep is in a 
gap, E z i should cause no net current along z' . Therefore 
the velocity averaged over on period driven by E z * is just 
the ratio of 8£ and the one period (= 5k z i /(e\E z > \/K)), 
which leads to v — (n/B)E z i x Sm. The Hall current 
due to E z i is calculated as 



G 



J = -pev 



—5m x E z >, 
ha 



(20) 



where p = eB/(ahn) is the density of states per subband 
and per unit volume. From this expression we can write 
the Hall tensor as cr±_ = — (e 2 /ha)5m±, where <Tjj = 
SijkO'k, and _L represents the component normal to z'. 
On the other hand, the normal component of the electric 
field, E± causes a classical drift in the direction normal to 
z' with the velocity v = (E ± x B)/B 2 = (n/B)E ± x 5m, 
which immediately leads to av = — {e 2 /ha)5m z >. Com- 
bining the two, we finally obtain the quantized Hall con- 
ductivity carried by the corresponding subband, 



<t = —(e 2 /ha)Sm. 



(21) 



We have thus derived the quantization of the 3D Hall con- 
ductivity from the fc-space hopping, as an approach alter- 
native to the usual Kubo formula. The problem of the 
transport in adiabatically varying potentials (so-called 
Thouless pumping) was first considered by Thouless 14 
for one-dimensional case. He showed that the 2D QHE 
in periodic potentials may be understood in terms of a 
fictitious ID system having two periodic potentials that 
slide adiabatically with each other. The discussion we 
have given here to describe the 3D QHE may be regarded 
as a two-dimensional version of the Thouless pumping. 

We can actually identify the Hall integers 5rm for each 
subband by keeping track of the motion of the orbits 
with k z i. Figure 2(a) typically depicts how a state in the 
first band moves by (1,0,0), by which we mean the state 
jumps once across a reciprocal cell boundary normal to 
k x (whose intersection is shown as one of the white lines 
in Fig. 2) but not across k y or k z . The state in the second 
band moves by (—1, 1,0). These triple numbers are the 
very Hall integers, (5m x ,5m y , Sm z ), carried by each sub- 
band (not to be confused with the total Hall integer) , and 
are in accordance with the small $ region in Fig. 1(a). 

We can now comment on the relation with the usual 
treatment of 3D QHE 3-5 . By summing up 5m ■ n = 1 
over the occupied subbands, we have m ■ n — r, where r 
is the number of occupied subbands and m is the sum- 
mation of dm over them. Since r is related to the filling 
of the tight-binding band in vb = r x p/(l/a 3 ), we obtain 



m. 



+ m v 6 v + m z 



(22) 



which coincides with the Diophantine equation (5) for 
s = 0. This means that in our picture above, only the 
gaps with s = are taken into account while the rest 
are neglected. We can actually show that the gap asso- 
ciated with s ^ only occurs, when the magnetic field 
is weak ($ <C 1), as a small gap between very disper- 
sive bands (versus k z i), so that the gap vanishes in the 
spectrum. Weak field regime is exactly the situation con- 
sidered here, since we are approaching to the quantum 
regime from the semiclassical one. The discussion here 
is for the lower half of the tight-binding band, while we 
can make a similar argument for the upper half in terms 



of hole orbits, which can be obtained as an electron-hole 
transformation (on the z/b-$ diagram) , which in turn cor- 
responds to s = — > 1 in the Diophantine equation. 

The mapped picture also gives an intuitive explana- 
tion why we have so few gaps for the symmetric case 
Fig.2(b),(c), and zero-component case (d). When two 
components in n coincide as in Fig. 2(b) with n oc 
(1, 1, 2), while two levels cross (for the value of k z * labeled 
as B), a gap does not arise because the couplings along k x 
and k v occur symmetrically (in a zigzag fashion) , so that 
the bands do not split, which is exactly the band touch- 
ing discussed in Sec. III. When the symmetry is even 
higher with n oc (1, 1, 1) in Fig. 2(c), couplings along k x , 
k y and k z all become symmetric and the band touch- 
ing occurs at every energy crossing. If n contains a zero 
component as in (d), two of the plane waves for s become 
parallel and the wells having the same depth become con- 
nected along the perpendicular direction. This results in 
a strong mixing between the states along that trough 
direction, so that the minibands for each value of k z > be- 
come wider. So, while the energy gap does arise when 
the energies of the adjacent troughs coincide, these gaps 
tend to be overlapped in energy by other, wide bands. 



(a) 1/9(1,2, 3) 




fes 2 rc/a 
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FIG. 2. Energy spectra versus k z i for Eq.(15) (left pan- 
els), where the final band structure is indicated as bars at 
the right-hand edge of each spectrum. Right panels are 
grayscale plots of the fictitious potential e(K x i , K y i, k z i) ver- 
sus (K x i , K y i ) (with lighter areas corresponding to dips) 
for different values of k z i (labeled as A, B, C in the left 
panels). Three sets of white lines represent the intersec- 
tions of the boundaries of the reciprocal unit cells normal 
to k x ,k y ,k z , respectively. We show three typical cases of 
(4>x,<t> v ,4>z) = (l/9)(n x ,n y ,n z ): (a) (n x ,n y ,n z ) = (1,2,3), 
(b) (1,1,2), (c) (1,1,1), and (d) (0,1,2). The energy 
spectra are periodic in k z i with period 2n/ (an) with 
n — \J n x + n y + n z . Motion of the orbits with k z i for the 
lowest subband (as well as for the second subband in (a)) are 
indicated, where the jumps to adjacent orbits are represented 
as arrows in the right panels and the resulting energy gaps 
as shaded ovals in the left. In (a) the energy spectrum for a 
weak field ((1/45)(1, 2,3)) are also shown for comparison. 

In the previous paper 12 we have shown that an ana- 
logue of Hofstadter's butterfly arises (i) as a function of 
the tilting angle 9 in magnetic field tilted in j/z-plane (ii) 
in anisotropic (t x » t y ,t z ) 3D lattices. If we apply the 
above argument to this case, we see that the fictitious 
potential e is dominated by the cosine band associated 
with the conductive direction, x. As shown in Fig. 3 this 
gives rise to one-dimensional bound states in its troughs 
and ridges with a small perturbation due to t y and t z . So 
the system can be regarded as an array of independent 
chains with two periods, which is described by a Harper's 
equation for the 2D Hofstadter's problem, and thus we 
can understand why the anisotropic case can have many 
gaps around the top and bottom of the band as seen in 
Fig. 3, while one component a yz is fixed to unlike the 
isotropic case. 
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FIG. 3. A plot similar to Fig. 2 for an anisotropic 3D sys- 
tem with 
(t x ,t y ,t z ) = (1.0,0.1,0.1) for {ct> x AvAz) = (1/15)(0,1,4). 
The contours highlight the troughs. 



VI. DUAL HALL CONDUCTIVITY 

We have just shown that two different systems related 
by the duality have similar energy spectra. It would be 
then interesting to see how the Hall conductivities are 
related between them for the corresponding energy gaps. 
For the 2D case, Thouless et aP have shown the rela- 
tionship between the Hall integers in the tight-binding 
and in the strong-field limits. Here we first review this 
in the context of our duality before moving on to the 3D 
case. Diophantine's equation for an energy gap in the 2D 
tight-binding case (a 2D version of Eq.5) is 



v B = s + t<j>, 



(23) 



where s,t are integers, vb the tight-binding band filling, 
and <f> the number of fluxes penetrating the 2D unit cell. 
The equation for the corresponding gap in the 2D strong- 
field case is written with the same s, t as 



v L = s + t{\/4>), 



(24) 



where vl is the Landau- level filling factor, i.e., the frac- 
tion of the filled states in the lowest Landau level, and 
4> is replaced with 1/0 according to the duality (Eq.9). 
From the relation vl = we can translate Eq.(24) 

into 



v B = s<j) + t. 



(25) 



Now we can see that the two topological integers s and t 
are interchanged between the dual [(tight-binding (23) <-> 
strong-field (25)] cases. The Widom-Streda formula, 10 ' 11 



indeed dictates that the Hall conductivity is the gradient 
in UB-(j> diagram, which gives — {e 2 /h)t in the former case 
— (e 2 /h)s in the latter. 

This results can be readily extend to the 3D case. 
Namely, if we take an energy gap in the 3D system with 
(cr yz ,a zx ,a xy ) = -f^(m x ,m v ,m z ), the Hall conductiv- 
ity in the corresponding 2D system becomes —(e 2 /h)s, 
where rrii and s are related via Eq.(5). 



VII. EXPERIMENTAL FEASIBILITY 



Let us comment on the magnitude of the magnetic 
field required to observe the energy gaps in 3D systems. 
For that it is essential that the coupling between the 
semiclassical orbits in the different wells exists. For the 
isotropic crystals considered here, the area enclosed by 
the semiclassical orbit in fe-space should be of the order 
of the typical well area, 27T/7 2 ~ (27r/a) 2 , which is simply 
Ba 2 /(/>() ~ 1. The situation is the same as in the 2D Hof- 
stadter problem, so the required field for atomic lattice 
constants is huge (B ~ 10 5 T for a =2 A). If we con- 
sider systems with larger unit cells, as in solid fullcrcne 
or zeolites with a =10 A, the required B is reduced to 
10 3 T but still large. In the anisotropic case discussed in 
Ref. 12 , by contrast, the required is much less stringent. 
This is because the typical well area is {2-K/a)d with d 
being the typical valley width, so the condition relaxes 
to 27r/Z 2 ~ (2ir/a)d, which can be made as small as one 
wishes by increasing the anisotropy (since t y , t z — > leads 
to d — > 0), although the scale of the energy gap shrinks. 
So it is a trade-off between large B required with large 
energy gaps (isotropic case) and smaller B suffices with 
small gaps (anisotropic). 

Another point is that a real 3D sample has always 
surfaces. Halperin and the present authors have shown 
in general that the 3D integer QHE should accompany 
quantized wrapping current, whose intensity and direc- 
tion are dictated by the quantum Hall integers 15 . This 
should apply to the present case of isotropic QHE. 



VIII. CONCLUSION 



We have investigated the energy spectra in the 
isotropic 3D lattice in magnetic fields applied in arbitrary 
directions, for which we have shown that the energy gaps 
arise unless the magnetic field points to high-symmetry 
directions. We have also calculated the quantum Hall in- 
tegers. In the latter part of the paper we have introduced 
a duality that relates B <-> 1/B (i.e., guiding center pic- 
ture in strong B and semiclassical picture in weak B) for 
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the present 3D problem, which gives the graphic expla- 
nation of the quantization of the Hall conductivity in 3D 
and the condition for the butterfly spectrum. 

M.K. would like to thank the JSPS for a financial sup- 
port. 
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